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Abstract 

The tunnel effect is considered liere witfiin the framework of elec- 
tromagnetic propagation. The classical problem of a plane gap of 
dielectric, surrounded on both sides by a medium with larger refrac- 
tion index, is studied in the case in which an electromagnetic plane 
wave impinges into the gap with an incidence angle larger than the 
critical angle. In this condition (total reflection), the gap acts as a 
classically forbidden region and behaves like a tunnel. The field inside 
the forbidden gap consists of two evanescent waves, each one having 
its wavefronts normal to the interface. In the present paper we study 
the total field derived as a superposition of two such evanescent waves, 
its wavefronts, and the directions of propagation of both phase and 
energy. 

In electromagnetism, an effect analogous to the tunnel effect of quan- 
tum mechanics occurs when a plane wave, propagating in a medium with 
refractive index n, impinges into a plane-parallel dielectric gap with refrac- 
tive index n' smaller than n, at an incidence angle larger than the critical 
angle io defined as sinzo = n/n'. If the thickness d of the gap is infinite, 
the field within the gap consists of a plane evanescent wave - attenuating in 
the direction normal to the interface - whose phase propagates in direction 
parallel to the interface (see Fig. la). If d is finite, the boundary conditions 
on both interfaces cannot be satisfied by a single evanescent wave and two 
evanescent waves, with the same direction of propagation of the phase, but 
attenuating into opposite directions, are required (see Fig. lb) [||]. 
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Figure 1: a) Border surface between medium 1, with refractive index n', and 
medium 2, with refractive index n < n' . For an incidence angle i larger than the 
critical angle iq = arcsin(n/n'), a single evanescent wave originates in medium 
2 and propagates parallel to the interface, b) Finite gap thickness of medium 
1, surrounded on both sides by medium 2. In this case, two evanescent waves, 
with the same direction of propagation but attenuating into opposite directions, 
originate inside the gap. The width of the gap d is shown, together with the 
coordinate system adopted in the theoretical analysis. 
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The properties of the total field inside the forbidden region are due to 
the fact that the two evanescent waves to be added have real amplitudes 
with opposite trends of variation in addition to the fact that the wavefronts, 
parallel to one another and normal to the interface, are not coinciding (on 
one wavefront the phase of a wave is different from the phase of the other 
wave, see Eqs.(||)). Accordingly, the phase of the total field varies along the 
wavefronts of the single waves to be added, and the wavefronts of the total 
field are not parallel to those of the component waves. 

Let us consider a system of Cartesian coordinates x,y,z (unit vectors of 
the axes i, j, k) with origin at O (Fig. lb), and an impinging TE wave with 
the electric field parallel to j, with direction of propagation s* = (ai +7k) in 
the plane xz. The ^/-component of the incident electric field can be written 
as 

= Eq exp[ikQn{ax + •yz)] , (1) 

where ko = u/c is the free-space wavenumber, Eq (which we assume to be 
real) denotes the amplitude of the incident field at the origin O and 

7 = Vl — q;2 . 

Inside the gap the total field is the superposition of two TE evanescent waves, 
whose electric field (still parallel to j) can be written as 

E~^ = pEq exp[ko{inax — Tz)] exp(— zcjt) , 

E~ = rEoexp[ko{inax + Tz)]exp{—iut) , (2) 

where 

r = V n'^a'^ — 1 

is a real quantity if, as assumed, the incidence angle is larger than the critical 
angle (a > l/n). 

The complex coefficients p and r can be deduced from the boundary 
conditions on the interfaces at z = and z = 
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From Eq. (^) it turns out that and ipr differ by 2$. All the above 
quantities are independent of the coordinates, and r represents the amplitude 
transmission coefficient of the gap[§]. The total electric field inside the gap 
can therefore be written as (apart from the time dependence exp(— icut)) 



E'^ = E+ + E- 



EoA\r\ 
2T 



exp{ikQnax) x 



exp[/cor((i — z)] exp{i(pp) + exp[—kQT{d — z)] exp{i(pj^ I 
and, by taking into account Eq.(^), we have 



E<^ 



EoA\t\ 
2T 



exp{ikonax + iifr) x 



exp[/cor((i — z)] exp(— z$) + exp[—koT{d — z)] exp(z$) 



(7) 
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From Eq.(l^) it turns out that the amphtude of depends only on z (the 
equi-amphtude surfaces are the planes z = constant), and decreases from the 



value £'oA|r|ycosh^(/i;or(i) — sin^ $/r, at 2; = 0, to the value Eo\t\, at z = d. 
The spatial dependence of the amplitude implies that geometrical optics is 
inadequate for describing electromagnetic propagation within a tunneling 
region 0. 

As to the phase, from Eq. it follows that the equation of wavefronts 
(equi-phase surfaces) is given by 

V9(x, z) = konax + \E'(z) + ip^- = cpo , (8) 
where <^o is a constant and 



tan[\&(2;)] = - tanh[/cor((i - z)] tan($) = - ^ tanh[kor{d - z)] 



cfj- = arctan 



2n-fT 



ta.nh.{kQrd) 



(9) 



Looking at Eq. (|), we can see that different wavefronts, corresponding to 
different values of the phase, are simply shifted in the x-direction with respect 
to one another. Figure 2 shows the wavefronts ipo = 0, vr, 27r for parameter 
values referring to an experiment dealing with frustrated total reflection in 
the range of microwaves 0. 

We can now derive both the wavelength and the equation of the rays. To 
this end, we have to evaluate gTa.d[ip{x, z)] and we have 

27r 

\giad[Lp{x,z)\\ 

and, for the ray equation x = X{z) defined as the lines of flux of gradyj, 

dX (dcp/dx) 



dz {dip/dz) 



11) 



From Eq. we have 



^ = kona (12) 
ox 
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Figure 2: Wavefronts (solid lines) as derived from Eq. (g), for tpo = 0, vr, 2tt, 
and rays (dashed lines) as derived from Eq. (|l^ , for three arbitrary values of the 
constant Xq. The parameter values refer to an experimental situation |Q] in which 
medium 1 consists of paraffin (n' = 1.49) and medium 2 of air (n = 1). Other 
parameter values are: a = 0.68, z/ = 10 GHz {to = liru), c = 30 cm/ns. 
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(15) 
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From Eq. (|T^) it follows that |grad<y9| > kgna and, from Eq. (jUlD, 

A < — < Ao, (na > 1) (16) 

na 

where Aq = 27r//co is the free-space wavelength (medium 2 in Fig. 1). We can 
conclude, therefore, that the total field inside the tunneling region is slow, 
that is the phase velocity along the rays is slower than the light velocity c. 
By integrating Eq. (p!5D, we obtain the ray equation 



X(z) 



z\-- 



A2 



4A;nr 



sinh[2fcor(ci - z)\ 



(17) 



where Xq is a constant. From Eq. (p!?]), we can see that (as was to be 
expected due to the symmetries of the problem) the rays are shifted in the 
x-direction, and that, since the amplitude of the total field is not constant 
with respect to 2, they are not straight lines (see Fig. 2). 

By denoting the angle between a ray and the 2;-axis with x (tan^ = 
dX/dz (see Eq. (|15|)), we have, at z = 0, 



smx 



A^ cosh^ {k^Vd) 



A2 cosh2(A;orci) - n272] + (^^l!) 



while, at z = d, 



sinx = a 

Since x is the refraction angle at z = and the incidence angle into the second 
interface a.t z = d, the unexpected conclusion is thus that the refraction low 
seems not to be valid for the rays inside the gap. 

The Poynting vector inside the gap has a component normal to the in- 
terfaces, contrarily to what happens for the field on the right of the single 
interface of Fig. 1. In order to evaluate the flux lines of the Poynting vector 
and their relationship with the flux lines of the phase, let us consider the 
vector S describing the energy propogation0 

S = iRe(EAH*), (18) 
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where the asterix indicates a complex conjugate. We easily obtain 



2Zo 
1 



-na 



'EoA\r\ 



cosh^[A;or(c/ - z)] - sin^ $ 



2Z, 



(19) 







where Zq is the free space impedance. As expected, the 2;-component of the 
Poynting vector does not depend on z. The flux lines of the Poynting vector 
can therefore be written as 



dx 
dz 



cosh.^[koT{d — z)] 



(20) 



and, by comparing this with Eq. (|T^), we see that the flux lines of the energy 
coincide with the flux lines of the phase. 

By means of Eq. (^, we are able to evaluate the phase difference Aip 
between two opposite points, O and O', along the z-direction inside the gap 
(see Fig. lb): we have 



A{p = {p{0, d) — (y9(0, 0) = arctan 



tanh(A;or(i) 



By including in the total phase also the temporal factor exp(— icut) disre- 
garded until now, the phase delay in going from O (at t = 0) to O' is 



= Aip/co. (21) 

In Fig. 3, we show as a function of the gap width, together with the time 
ti = d/c. If we suppose to perform an experiment (like the one reported in 
Ref. 0) with a monocromatic wave and we put two probes at O and O', we 
would actually measure a phase delay as given by Eq. (|21|) which is, without 
doubt, an observable To clarify the physical meaning of this delay we 
have to complete the analysis by considering not a plane wave impinging the 
gap, but a narrow beam or a wave packet]]^, ^. However, this improvement 
exceeds the purpose of the present work and will be presented elsewhere. 

Finally, we wish to note that becomes independent on d, for large 
d (see Eq. (pT])). A behaviour of this kind was also obtained within the 
framework of a quantum-mechanical theoretical model, due to HartmanQ, 
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Figure 3: Phase delay time along the z-direction as derived from Eq. (|2l|), 
together with the time ti = z/c, as a function of z. Parameter values are the same 
as those in Fig. 2. 

for a particle tunneling through a rectangular potential barrier. Also in that 
case, the "traversal time" under barrier tends to be constant for large barriers, 
and the superluminal effect so obtained is known as "Hartman effect". It is 
not easy to understand the nature of that time but its behaviour, for large 
barrier, very similar to the one as derived here (Eq. (|2l|) ), could characterised 
it as a phase-delay 1^, p!0| . 
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